Fractional quantum Hall state at v = 1/4 in a wide quantum well 
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We investigate, with the help of Monte-Carlo and exact-diagonalization calculations in the spher- 
ical geometry, several compressible and incompressible candidate wave functions for the recently 
observed quantum Hall state at the filling factor v = 1/4 in a wide quantum well. The quantum 
well is modeled as a two-component system by retaining its two lowest subbands. We make a di- 
rect connection with the phenomenological effective-bilayer model, which is commonly used in the 
description of a wide quantum well, and we compare our findings with the established results at 
v = 1/2 in the lowest Landau level. At v — 1/4, the overlap calculations for the Halperin (5,5,3) 
and (7,7,1) states, the generalized Haldane-Rezayi state and the Moore-Read Pfaffian, suggest that 
the incompressible state is likely to be realized in the interplay between the Halperin (5,5,3) state 
and the Moore-Read Pfaffian. Our numerics shows the latter to be very susceptible to changes in 
the interaction coefficients, thus indicating that the observed state is of multicomponent nature. 

PACS numbers: 73.43.Cd, 73.21. Fg, 71.10.Pm 



I. INTRODUCTION 

Advances in fabrication of high-quality GaAs semicon- 
ductor systems have led to an ever growing collection 
of the observed incompressible fractional quantum Hall 
states in a variety of settings* 1 - These states occur at par- 
ticular ratios between the number of electrons N and the 
number of magnetic flux quanta that pierce the sys- 
tem in the direction perpendicular to the sample. This 
commcnsurability can be expressed as the filling factor 
v = N/Nff, = p/q in terms of integers p, q, which is 
the single most important quantity that characterizes the 
quantum Hall state. 

In a thin layer, q usually turns out to be an odd in- 
teger, the fact which had its pioneering explanation in 
terms of the Laughlin wave function 2 - for the case of 
p = 1, q = 3, 5, 7, ... and its subsequent generalizations in 
terms of composite fermions 3 - (CF), applicable to general 
integers p, q as long as q is odd, and hierarchy theory* 4 - 
However, a state with an even denominator has also been 
observed^ but in the first excited Landau level (LL). One 
cannot account for it in the usual Laughlin/compositc 
fermion approach and the idea of pairing has commonly 
been invoked to explain the origin of this fraction^ 
The simplest realization of pairing between spin polar- 
ized electrons is the so-called Pfaffian defined by the 
Moore-Read wave function 7 - and supporting excitations 
with non-Abelian statistics^ 

The possibility of an extra degree of freedom lifts the 
requirement of Fermi antisymmetry and hence gives an- 
other route towards realizing even denominator fractions. 
The additional degree of freedom can be the ordinary spin 
or else a "pseudospin" in case of a wide quantum well, 
where the two lowest electronic subbands correspond to 



t, i- If the sample is etched in such a way to create a 
barrier in the middle, thus supressing tunneling between 
the two "sides" , one can think of it as a bilayer with | , J, 
denoting the left and right layer where electrons can be 
localized. Incompressible quantum Hall states for such 
systems have been theoretically predicted in Ref. and 
experimentally confirmed for cases of bilayer at filling 
factor v — 1 and v = 1/2 . 10 ' 11 Later on, essentially the 
same quantum Hall state at v = 1/2 was observed in a 
sample which had the geometry of a single wide well pi 2 - It 
was argued, on the basis of a self-consistent Hartree-Fock 
approximation, that in a wide well the electrons (due to 
their mutual repulsion) reorganize themselves so as to 
form an effective bilayer distribution of charge. Hence, 
an equivalence between the two very different samples 
was claimed and theoretical works set out to analyze the 
problem from this premise d 3 ' 14 

On the basis of the quantum mechanical overlap with 
the ground state obtained in exact diagonalization (ED), 
including a realistic bilayer confinement potential, Ll.ef.ll3l 
established that the ground state is well described by 
the so-called (3,3,1) Halperin wave function^ This wave 
function distinguishes between two kinds of electrons and 
the fact that it describes the system is what we mean by 
the system being "multicomponent" . Experimental work 
gave further insight into the nature of the multicompo- 
nent state at v = 1/2 and strengthened the belief that the 
(3,3,1) wave function is a correct physical description^ 
Namely, the behavior of the excitation gap as a function 
of tunneling amplitude As^s (i.e. the splitting between 
the two lowest subbands) was found to have upward cusp 
at the intermediate value of As^g and the state was 
quickly destroyed by the application of electrostatic bias 
(charge imbalance) ^ In Ref. [HI, a numerical study was 
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able to reproduce the observed upward cusp in the acti- 
vation gap by diagonalizing the bilayer Hamiltonian with 
explicit intcrlaycr tunneling. 

A recent experimental papei— reports the observation 
of the v = 1/4 quantum Hall state in a wide quantum 
well. The state is fragile and almost indiscernible when 
only a perpendicular magnetic field is applied (although 
one could expect that with yet higher sample qualities, a 
small plateau would be developed already at that point). 
However, when the magnetic field is tilted, there is a clear 
dip in the value of longitudinal resistance R xx , signifying 
the presence of an incompressible state. 

In this paper we analyze the complex interplay between 
the single- and multicomponcnt nature of the ground 
state at v = 1/4 in a wide quantum well, in comparison 
with the ground state at v = 1/2. Contrary to previ- 
ous studies ; 13 ! 14 we do not make the ad hoc assumption 
that the wide quantum well may be described as an ef- 
fective bilayer. Instead, we consider the two lowest elec- 
tronic subbands of the quantum well, which is modeled 
by the infinite square well for the sake of convenience 
but cross-checked with other confinement models. The 
energy splitting between these two subbands, the asso- 
ciated wave functions of which are symmetric and anti- 
symmetric, respectively, in the z-direction, is given by 
Asas (occasionally referred to as the tunneling ampli- 
tude). Due to the low filling factor (y = 1/4), the power 
of the ED method will be rather limited and other com- 
plementary approaches may be needed to fully explain 
the experimental findings. 

The paper is organized as follows. Sec. |TT] is devoted 
to the single-component candidate for v = 1/4, and we 
study its overlap with the exact Coulomb ground state 
within various confinement models. In Sec. IIII1 we de- 
fine the multicomponcnt wave functions expected to be 
relevant at this filling factor. The two likely candidates, 
the Halpcrin (5,5,3) and (7,7,1) states, are investigated 
within a simple bilayer model without tunneling. The 
two-subband model of the quantum well is introduced 
and described in Sec. IIVI Our main results of ED calcu- 
lations in the spherical geometry are presented in Sec.lVl 
To extend the reach of our numerics, we furthermore de- 
ploy Monte-Carlo simulations of the trial wavefunctions 
identified beforehand to analyze their energetic competi- 
tion. We summarize with our view on the nature of the 
state at v = 1/4 in Sec. [VTJ 



II. ONE COMPONENT STATE 

A. Pfafflan at v = 1/4 

There is a natural candidate for the fully polarized 
quantum Hall state at v = 1/4 - it is the generalized 
Moorc-Rcad Pfaffiani^ 



#P f (zi,...,.zjv) = Pf 



expressed in terms of the complex coordinate of the elec- 
tron in the plane where zj = Xj + iyj. The object Pf is 
defined as 
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acting upon the antisymmetric N x N matrix My, and 
Sn is a group of permutations of N objects. Pf renders 
the wave function totally antisymmetric and encodes the 
same kind of correlations as in the more familiar v = 5/2 
c&se£ In the spherical geometry 4 ^ many-body states are 
characterized by the number of electrons N, the number 
of flux quanta generated by a magnetic monopolc 
placed in the center of the sphere and extending radially 
through its surface, and an additional topological number 
which is the shift. For the Pfaffian in Eq. (fTJ), the three 
numbers are related by the formula = 47V— 5. \I/pf is a 
zero-energy eigenstate of a certain 3-body Hamiltonian,- 
but in our calculations it was generated from its root con- 
figuration via the squeezing technique^ On the other 
hand, the Coulomb (two-body) Hamiltonian commutes 
with the angular momentum operator L because of ro- 
tational invariance and, by Wigncr-Eckart theorem, the 
interaction is parametrized by discrete set of numbers Vl 
known as the Haldane pseudopotentials^ The motion of 
electrons is therefore fully described in terms of the in- 
plane (spherical) coordinates 8, 4> and the use of different 
confinement models in the (perpendicular) z-dircction 
(neglecting the in-plane magnetic field) will only mod- 
ify the values of pscudopotcntials. 



B. Finite thickness models 

Most of the candidate wave functions for quantum 
Hall fractions have been extensively studied via numeri- 
cal techniques such as ED or Monte Carlo. For the sake 
of convenience, but also due to the intrinsic ambiguity 
which stems from the fact that in a strongly correlated 
system many input parameters (e.g. the precise form of 
the interaction) are unknown, it is natural to start off 
from the limit of infinitely thin layer of electrons inter- 
acting via Coulomb force and hope that the inclusion 
of e.g. realistic confinement and sample thickness will 
have small, perturbative corrections. There have been 
different proposals to account for the finite thickness of 
the sample in the perpendicular direction, but the one 
that is straightforward and most natural from the point 
of view of ED is the Zhang-Das Sarma (ZDS) models 
which is simply given by substituting the interaction 
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(we will always denote by w the width of the sample and 
the energy is always expressed in units of e 2 /elg, where 
the magnetic length is Is = \JhcjeB is given in terms of 
the perpendicular magnetic field B). Qualitatively, this 
substitution softens the interaction^ and was studied ex- 
tensively (together with other confinement models, some 
of which we will introduce below) in Ref. HO, where it 
was advertised to significantly stabilize the Moore-Read 
Pfaffian at v — 1/2 (the effect being most pronounced in 
the second LL), but (in most cases) decrease the over- 
lap somewhat for the Laughlin states at v = 1/3 and 
1/5. In Ref. [2l| it was noticed that this kind of interac- 
tion can lead to an instability of the composite fermion 
sea, which is believed to describe the compressible state 
at v = 1/2 in the lowest LL, towards the paired state 
described by the Pfaffian. Indeed, the CF Fermi liquid 
can be regarded as a special member of the general class 
of paired CF wavefunctions^ of which it represents the 
limit of vanishing gap. 

Although the ZDS model © has a very simple form, 
there is no physical wave function that corresponds to 
this confinement potential in the z-direction. Other pop- 
ular choices for the confinement in the ^-direction include 
the infinite square well (ISQW) and Fang-Howard (FH), 
which are presumably more realistic than ZDS because 
they arc defined by the actual wave functions of simple 
model potentials for the quantum well, given by 
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respectively 



0FH(Z) = VP 



C. Overlaps 



-3z/2u 



(3) 
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We have performed ED calculations for various con- 
finement models and all system sizes N = 
6,8,10,12 accessible at present. In Fig. [1] we present 
the overlap | (^pt Inexact) | between the exact Coulomb 
ground state at v = 1/4 and ^pf, finite width being 
modeled by the ZDS ansatz ©. The size of the Hilbcrt 
space at N = 12 is noteworthy: the dimension of the 
L z = sector is 218 635 791. 

It appears that the overlap of the Pfaffian state is 
rather high for large values of the width (even if it is 
negligible for small w's). These values could likely be 
increased further by considering general pairing wave 
functions^ 2 - However, these overlaps alone cannot be 
taken as solid evidence for a pairing nature of the v = 1/4 
for two reasons. First, for N = 6 and TV = 12 there is 
the aliasing problem with composite fermion states: Jain 
states with different physical properties (e.g. Abclian 
instead of non-Abelian statistics) occur at the same val- 
ues of N and on the sphere (because of finite system 
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FIG. 1: (color online) Overlap |(^/pf Inexact )| between the ex- 
act Coulomb state for finite width (ZDS model) and the Pfaf- 
fian at v — 1/4. 



size). High overlap for the aliased states may therefore 
come from other incompressible states different from the 
Pfaffian. Secondly, for the non-aliased states at N = 8 
and N = 10, there appears to be a critical value of the 
width at which the overlap as a function of w suffers a 
sharp jump. By analyzing the entire low energy spectrum 
on the sphere as a function of width, we have established 
that the (neutral) gap collapses at the critical point of 
w/Ib- Therefore, in order to get to the Pfaffian phase, 
one must go through a (first order) phase transition. Be- 
fore the transition, the ground state is obtained in the 
L > sector of the Hilbert space and the overlap with 
the Pfaffian (which resides in L = sector) remains zero 
due to the difference in symmetry. 

The lack of adiabatic continuity and the aliasing prob- 
lem cast some doubt on the Pfaffian state as a good can- 
didate for v = 1/4 in the lowest LL. We have also checked 
using other confinement models ([21 E]) , but in these cases 
for N = 8 and N = 10 the overlap remains zero for any 
value of w/Ib ■ Thus our ED results do not yield a def- 
inite answer with respect to the relevance of "fpf in the 
single layer at v = 1/4. 

We would like to stress the qualitative difference in 
our results obtained by using ZDS versus other confine- 
ment models which appears, to the best of our knowl- 
edge, to be the first such case in the literature. The 
smaller overall energy scale (and the smaller gap as well) 
is very likely to be at the origin of this discrepancy. We 
note in passing that, contrary to the finite- width mod- 
els which change all pseudopotentials at once, one may 
start from the pure Coulomb interaction and vary just a 
few strongest pseudopotentials.— We have tried varying 
both V\ and V3 , but this procedure does not stabilize 
the Pfaffian phase in any finite region of the parameter 
space for N = 8. 
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III. TWO-COMPONENT STATES 



tonian 



Soon after Laughlin's wave function describing the in- 
compressible state at v = 1/3 when the electron spins are 
fully polarized, Halpcrini^ proposed a class of generalized 
wave functions defined as 



^mm'n (^1 ) •••) ! ^1 ' "*' ^ N ± ) 

= n^ T - z )r - z ir' niN - z *) n > ( 5 ) 

i<j k<l s t 

where the electrons are distributed over two compo- 
nents (labeled by Til)- The exponents m, ml denote the 
"intra" -component correlations originating from the ba- 
sic Laughlin-Jastrow building blocks within each com- 
ponent, whereas n describes "inter" -component correla- 
tions (we have omitted the ubiquitous Gaussian factors 
and implicitly assume that there is a spinor part to this 
wave function as well as an overall antisymmetrization 
between f and [). In order for these wave functions to 
be eligible candidates for the ground state of the system, 
one must enforce an additional requirement that they be 
eigenstates of the Casimir operator of the SU(2) group, 
i.e. the total spin S 2 , as long as the interaction is sym- 
metric with respect to "intra" and "inter" components 
(e.g. the usual case of electrons with spin). However, 
apart from electrons with spin, the wave functions (J5]) 
have also been used in bilayer systems where this sym- 
metry is broken as soon as the layer separation is non- 
zero. In this case, the wave functions ([5]) need not be 
eigenstates of the total spin. There have been general- 
izations of these wave functions in the physics of bilayer 
systems at total filling facto r 24 i 25 i 26 v = 1 and to more 
than two components^ where further constraints on the 
possible values of to, to', n were derived within the plasma 
analogy^ In a two-component case, these turn out to be 
the intuitive requirement that "intra" -component inter- 
actions are stronger than "inter" -component interactions: 
to, to' > n. For the particular case of two components 
and to = ml = n + 2 (which includes ^331 and ^553), the 
Halperin wave function ^ can be analytically cast into 
a paired for m 26 ' 29 via Cauchy determinant identity (up 
to the unimportant phase factor), 
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where the pairing function is given by det 



In the 



case of the Ill-state, this pairing nature was recently ex- 
ploited to make a connection to paired composite fermion 
states, and to construct wave functions interpolating be- 
tween these two regimes^ Halperin wave functions are 
the exact zero-energy eigenstates of the two-body Hamil- 



= E 



J2 VpPjJiN* -L)+J2 V^P^N* - L) 



L=0 

n-1 



L=0 
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where P^ a (L) projects onto the state with angular 
momentum L of particles i and j with respective 
(pseudo)spins a and a'. Besides offering great conve- 
nience for handling Halperin wave functions ((5]) in ED, 
Eq. ^ enabled counting of the number of excited quasi- 
hole states and reaffirming the idea that the states de- 
scribed by Eq. © possess Abelian statistics^ 

At the filling factor v = 1/4, there are three wave 
functions of the form ([5]) that meet the necessary phys- 
ical requirements, ^553 = (5,5,3), ^771 = (7,7,1) and 
^5131 = (5,13,1). None of them is an eigenstate of 
S 2 , so they are more adapted to the case of a bilayer 
than that of real spin. In Fig. [2] we present the ba- 
sic overlap characterization of the first two wave func- 
tions in a simple bilayer model defined by the interaction 
W T (r) = V ll (r) = l/r,V n (r) = 1/V'r 2 + <P (d being 
the distance between the layers)^ (5, 5, 3) displays a fa- 
miliar maximum in the overlap for small distance be- 
tween the layers. (7, 7, 1) was dismissed in Ref. [lfl argu- 
ing that it would more likely lead to two coupled Wigner 
crystals than an incompressible liquid. Our diagonaliza- 
tion scheme is not adapted to address states with broken 
translation symmetry, so we do not see an a priori reason 
to reject this state. The results in Fig. [2] are for N = 8 
particles, they are fully consistent with those of smaller 
N, but direct comparison between (5,5,3) and (7,7,1) 
is not possible because they are characterized by differ- 
ent shifts (—5 and —7 respectively)^ We will address 
this issue below by extrapolating to the thermodynamic 
limit the respective trial energies from Monte-Carlo sim- 
ulations for both of these states. 

The last possibility, (5, 13, 1), is a peculiar one because 
it can only occur in the case of a strong density imbal- 
ance. Such an imbalance would lead to an increase in 
the charging energy, but if one of the coupled states is a 
prominent quantum Hall state, the gain in correlation en- 
ergy can outweigh the price of charge imbalance, as it has 
been experimentally verified^ However, in the present 
case, our numerical calculations confirmed that this can- 
didate can be discarded because it takes unrealistically 
high values of the sample width for this wave function to 
have any numerical relevance at all. 

Given the low filling factor v — 1/4 we are studying, 
one must also consider the possibility of nearby compress- 
ible states that can intervene for some values of the ex- 
ternal parameters. Apart from the obvious metallic state 
similar to the Fcrmi-liquid-likc state proposed by Rczayi 
and Read)22 there is in principle also the Haldane-Rczayi 
(HR) state which is defined by 
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FIG. 2: (color online) Overlap between the exact bilayer state 
with (5,5,3) and (7,7,1) states for N = 8 particles. 



*HR(K T ,^}) = dct 
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The last term is a global Laughlin-Jastrow factor for 
all particles regardless of their spin. \I/hr. is the zero en- 
ergy eigenstate of the interaction parametrized by the set 
of pseudopotentials Vl = {1,1,0, 1,0,...} and occurs at 
the shift of —6. It is also a spin singled and compress- 
ible on the basis of its nonunitary parent conformal field 
theoryi£i22 However, its edge theory^ is closely related 
to that of the Abelian (5, 5, 3) state, which suggests that 
the HR state may be in the vicinity of the incompress- 
ible state and nonetheless affect the physical properties 
of the system. Recently there have been proposals that 
compressible states can be molded into incompressible 
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and the Vs represent monopole spherical harmonics with 
—N^/2 < m < N<p/2 (wc assume that the states are 
entirely within the lowest LL). We refer to the states 
and @ as symmetric and antisymmetric, respectively, 
because of their reflection symmetry with respect to the 
center of the well. If their energy difference is denoted by 
Asas, the corresponding second quantized Hamiltonian 
is given by2£ 
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(9) 



rn. 2 , r?i3 , 7714 1 cr 1 ^m^a^ ^m^a^ ^m^a^ ; 



{m} {a} 



where cill annihilates (creates) an electron in the state 
m with pseudospin a. 

The matrix elements V^'^^fm, can be straightfor- 
wardly evaluated from the Haldane pseudopotentials for 
the resulting in-plane interaction 
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IV. THE QUANTUM- WELL MODEL 

So far we have discussed the stability of the one- 
component Pfaffian state in different finite-width mod- 
els (Section |n|, and two-component states in a bilayer 
model where each layer is considered as an infinitely nar- 
row quantum well (Section IIIip . In this Section, we con- 
sider an infinite square well of width w in the direction 
z € [0,w]. The electronic motion in the z-direction will 
then be quantized, yielding an electronic subband struc- 
ture. 



A. Two-subband approximation 

Instead of a full description with all the electronic sub- 
bands, we only consider the two lowest subbands and 
identify them with the two pseudospin states, ^f.j. = 
( l>hl( z ) Y N 4 ,/2,N 4 ,/2jn( ' ( t>)' where 



where the position variables are expressed in units of Is 
such that the integral is dimensionless. 

In this paper we do not make an attempt to quan- 
titatively model the experiment of Ref. [T5 , but we are 
interested in the possible phases that may occur and 
the transitions between them. Therefore, we expect the 
model described by the Hamiltonian © to be qualita- 
tively correct and in agreement with other confinement 
models that assume the lowest subband to be symmet- 
ric and the first excited one to have a node in the centre 
(z = w/2). Any difference of the confining potential away 
from the infinite square well will modify the energy eigen- 
values and the associated wave functions 4>&(z). However, 
it is expected that the energies arc more strongly affected 
than the wave functions. In particular, the nodal struc- 
ture of the wave functions is robust, such that the two 
lowest eigenstates of the infinite well faithfully represent 
the underlying features. However, we will allow for the 
general values of the level splitting Asas to account for 
the variations in the eigen-energics. 
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B. Connection between the quantum- well model 
and the bilayer Hamiltonian 

From a more general point of view, the quantum-well 
model exposed above is a two-component model such as 
the bilayer model, which has been used in the discussion 
of the wide quantum well.— Indeed, the wide quantum 
well allows the electrons to reduce their mutual Coulomb 
repulsion by exploring more efficiently the z-direction, 
and it has been argued that due to this effect, a sponta- 
neous bilayer may be formed, under appropriate condi- 
tions, in a wide quantum welL 12 ' 13 Here, a connection is 
made between both two-component models, on the basis 
of the Hamiltonian ©. The intermediate steps in the 
derivation of the effective model may be found in the 
Appendix A. 

The Hamiltonian © may be rewritten in terms of the 
density and spin-density operators projected to a single 
Landau level. The Fourier components of the projected 
density operator of pseudospin-cr electrons reads 



'(<0 = E ( 



-zq-R. I 



in terms of the 2D wave vector q and the guiding-center 
operator R, the latter acting on the states labeled by the 
quantum numbers m. It is furthermore useful to define 
the total (projected) density operator 

p(q) = p T (q) + ^(q) (ii) 

and the projected pscudospin density operators, 



5"(q) 



E 

m.tn' 



(m \e lqR | r, 



A t a ' a 



(12) 



where r M , are the usual 2x2 Pauli matrices with a = 

<7,<7 r 

x,y,z. 

In terms of the projected (pseudospin) density opera- 
tors, the Hamiltonian §§§ approximately reads 



q 

+2^^(q)^(-q)^(q) 



-AsasS z (cl=0) , 



(13) 



where the SU(2)-symmetric interaction potential 
Ksc/(2)(q) and the symmetry-breaking potential V^(q) 
are linear combinations of the Fourier-transformed 
potentials defined in Eq. (|10[) . Their precise form 
is given in the appendix by Eqs. (|A2[) and (|A5|) . 
respectively. The Hamiltonian (|13|1 neglects a particular 
term oc S z (— q)S z (q), which turns out to constitute the 
lowest energy scale in the interaction Hamiltonian ([9]) 
[see Eq. (|A9|) in the appendix]. 
Furthermore, 
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cIb Ib 



(14) 



is the effective subband gap. The numerical prefactor 
7 depends on the precise nature of the considered con- 
finement potential, and, as shown in the appendix, the 
expression (|14[) is derived within a mean-field approxi- 
mation of a particular term in the Hamiltonian ([9]) . The 
expression (|14j) is easy to understand - whereas the sub- 
band gap Asas tends to polarize the system in the | 
state, namely in narrow samples, the second term in 
Eq. (p~4|) indicates that the interactions are weaker in 
the | subband. From the interaction point of view, it 
is therefore energetically favorable to populate the first 
excited subband. This effect becomes more pronounced 
in larger quantum wells. Notice furthermore that this 
argument also delimits the regime of validity of the two- 
subband approximation of the wide quantum well; when 
the term jv(e 2 /cIb) x (w/Ib) becomes much larger than 
the bare subband gap Ag^g, the electrons may even pop- 
ulate higher subbands, which are neglected in the present 
model, and the system eventually crosses over into a 3D 
regime. 

Notice that the Hamiltonian (fT5|) has the same form as 
the Hamiltonian which describes a bilayer quantum Hall 
system^ up to a rotation from the z- to the £-axis. In 
this rotated reference frame, one may define the intra- 
and inter-layer interactions as 



V A (q) = V su{2} (q) + V* b (q) 



(15) 



V, 



-V, 



2 T D TTT (q) 
2 T D TU (q) 



and 



V E (q) 



Vsu(2) (q 
l 

V 2D 



V£(q) 



(16) 



^ TTT (q) 



UU, 



-v. 



2D 

As for the case of the true bilayer, the thus defined intra- 
layer interaction is stronger than the inter-layer interac- 
tion, for all values of q. 

Since our ED calculations employ the Hamiltonian 
in order to compare the numerical results with the 
Halpcrin states ([5]) which are the native eigenstates of 
the true bilayer Hamiltonian ©, we can apply the map- 
ping between the two models described above in a re- 
verse fashion. As Halperin wave functions are commonly 
labeled by the single particle states | f), I I) (which are 
the eigenstates of S z ) and defined by interaction poten- 
tials {Va,Ve}, we can imagine a linear transformation 
(rotation from z to x) that transforms them into (unnor- 
malized) symmetric |+) = | T) + I I) an< l antisymmetric 
|— ) = | t) — | J.) combinations. Then, by inverting the 
Eqs. (fT"5)) and (fll)]) , we obtain the set of interaction po- 
tentials that generate the Halperin states (m,m',n) in 
a quantum well description. In what follows, Halperin 
states (O are understood to be indexed by |+), |— ) in- 
stead of the usual notation | |), | I), unless explicitly 
stated otherwise. 



7 



C. Energetics of trial wavefunctions 

To extend the reach of our calculations to system sizes 
larger than those which can be treated in ED, we set up 
Monte-Carlo simulations of the trial states which have 
emerged as good candidates for the ground state. The 
general strategy of this approach is to obtain an estimate 
of the energy in the thermodynamic limit for the different 
trial states based on a scaling with system size of their 
energies. 

As detailed in section IIVBI above, we expect forma- 
tion of two-component wave functions where S x is a good 
quantum number, such that the Halperin wave functions 
are expressed in terms of the coordinates of electrons in 
the |+) and |— ) states, and lower well, indexed below 
by a. We consider cases with equal population of elec- 
trons in these two bands, or full population of the lowest 
subband in the ISQW for the single component cases. 

In order to calculate efficiently the interaction of elec- 
trons in a well of finite width using Monte-Carlo simu- 
lations, we replace the interaction ([10)) with an effective 
potential that reproduces all pseudopotential coefficients 
of the original potential V2D- Many such potentials can 
be constructed. Here, we use an interaction of the form 
proposed in Ref. l37l . built from simple polynomials^ 

V:/{r)=Y J cfr k . (17) 
k=-l 

The pseudopotentials of the monomials r n can be evalu- 
ated analytically (generalizing Ref. [H) . Choosing Ck to 
match the pseudopotential coefficients of the interaction 
(|10p becomes a simple linear problem. Crucially, we allow 
for the coefficient of the Coulomb term c_i to be varied, 
also. The number of terms is chosen equal to the minimal 
number required to match the relevant pseudopotentials 
(odd pseudopotentials V^m+i for intra (pseudo-)spin in- 
teractions and all N$ + 1 terms, otherwise). 

It is habitual in the literature to introduce a neutral- 
izing background, in order to highlight the correlation 
energy associated with a wavefunction. We use a back- 
ground Eh g [4>] that matches the distribution |0(z)| 2 of 
electrons in their subbands, in order to study the correla- 
tion energy of the different states. However, to establish 
a final comparison between the different wavefunctions, 
a unique convention for the background is required, and 
we adopt the background of the single layer configuration 
as a reference point E re { = E\, s [<fif]. 

Extrapolation to the thermodynamic limit is under- 
taken as two separate steps. The correlation energy is 
obtained by linear scaling over the inverse system size 
N _1 , using the habitual rescaling of the magnetic length 
l' B = [iV/^A^)] 1 / 2 ^.— For the two-component states, 
the difference in background energy E Te { — Eg [^1] 
is extrapolated separately, and added to the correlation 
energy. 



V. COMPETING PHASES IN THE 
QUANTUM- WELL MODEL 

In order to justify the model of the quantum well, in 
this Section we present the ED study of the Hamiltonian 
© and analyze the energetics of the relevant trial wave- 
functions in Monte-Carlo simulations. We briefly revisit 
the problem of v = 1/2 extending the results of Refs. [l2l 
andll4l (Section |VA[) and then present results pertaining 
to v = 1/4 (Section EH). 

A. i/ = 1/2 in a quantum well 

At this filling factor the competing phases we consider 
here are the (3,3,1) Halperin state, the Moore- Read Pfaf- 
fian and the HR state. Ref. [U has demonstrated a com- 
petition between the multicomponent (3,3,1) state and 
the fully polarized single-component Moore-Read Pfaf- 
fian. In the region of small tunneling, the ground state 
shows high overlap with the Halperin state; as the tun- 
neling is increased, the Halperin state is destroyed and 
the Pfaffian takes over. The point of crossover between 
the two is related to the upward cusp in the activation 
gap-i^ 

Fig.[3]shows our ED results for 8 particles in the quan- 
tum well at the filling factor v = 1/2. Figs. OJa), (b), and 
(c) represent the overlap between the exact ground state 
and the (3,3,1), the Pfaffian, and the HR states, respec- 
tively, as a function of the well width w/Ib and the bare 
subband gap As as- In general, the latter is a mono- 
tonically decreasing function of the well width. Again, 
we choose w and As^s as independent parameters of 
the model. Furthermore, we plot the quantity denoted 
by (S z ), the expectation value of the S z = S z (q = 0) 
component of the pseudospin which has the meaning of 
the "order parameter" [Fig. |3fd)]. One notices that (S z ) 
continuously crosses over from a full polarization in the 
"T subband at low values of w/Ib and a large gap Asas 
to a polarization in the J, subband for larger quantum 
wells and small gaps As as- As it is discussed in the 
previous section, the interactions in a wider quantum 
well favor a population of the first excited electronic sub- 
band I because of the node in the wavefunction in the 
z-direction and, therefore, decrease the effective subband 
gap. Indeed, Eq. (fT4|) indicates that the crossover line 
from positive to negative As^s is characterized by a bor- 
der that is linear in w/Ib- This behavior is also apparent 
in Fig.[3fd). Notice, however, that for large negative po- 
larizations (large negative As^s), the two-subband ap- 
proximation is no longer valid and the occupation of even 
higher electronic subbands must be taken into account, 
as already mentioned in Sec. IIVBI 

Note, furthermore, that we have defined our (3,3,1) 
state to be an eigenstate of the S x operator in the termi- 
nology of the true bilayer and not the usual S z operator 
(naively defining the Halperin state to be the eigenstate 
of S z does not give any appreciable overlap with the exact 
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FIG. 3: (color online) Overlap between the exact Coulomb 
state of the quantum well for N = 8 particles at v = 1/2 
with the Halperin (3,3,1) state (a), the Pfaffian (b) and HR 
(c) states. The expectation value of the S z component of the 
pseudospin is plotted in (d). 



ground state). There is a simple reason why this needs 
to be done: because the states of the quantum well pos- 
sess nodal structure (0, the true bilayer states (like the 
Halperin states) need to be rotated first from z to x di- 
rection, in order to match this symmetric/antisymmetric 
property, before direct comparison can be made. 

With this convention, the (3,3,1) state has its largest 
overlap (< 0.95) with the exact ground state in the vicin- 
ity of the crossover line (S z ) = 0. However, the over- 
lap remains quite large even in regions beyond this line, 
where the polarization becomes non-zero [Fig. [3Ja)], in 
agreement with Ref. H3 . This behavior may have two 
different origins. First, one notices that S z is not a good 
quantum number if the SU(2) symmetry-breaking terms 
of the Hamiltonian (|13p in the x-direction are taken into 
account. Especially in the vicinity of the crossover line 
As^s ~ 0, the symmetry breaking is governed by these 
terms in the x-direction, and S x (q = 0), which does not 
commute with S z , is expected to be a good quantum 
number. An alternative origin of the large overlap with 
the (3,3,1) state even in regions with (S z ) ^ may be a 
possible admixture (~ 5%) of states to the ground state 
that are orthogonal the (3,3,1) and possess a finite po- 
larization in the z-dircction. 

The largest values of the overlap between the compress- 
ible HR state and the exact ground state are also found in 
the vicinity of the crossover line from positive to negative 
As as, though at extremely large values of w/Ib- Notice 
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FIG. 4: (color online) Energies in the thermodynamic limit 
for the (3,3,1) and Pfaffian states at v = 1/2 (left). Data 
shown is for the infinite square well as a function of the well 
width w. The correlation energies are shown with respect to 
the single component background. In the absence of tunnel- 
ing, the (3,3,l)-state has lower energy at all w. The critical 
tunneling strength required to favor the Pfaffian state (top 
right) and a few typical differences in the extrapolation of 
the background energies for different values of the well width 
(bottom right). 



that the overlap (0.64 for w/Ib — 10.0) is generally much 
lower than for the (3,3,1) state. At large values of the 
bare subband gap Ag^s (and narrow quantum wells), 
the system becomes polarized in the t subband, and the 
ground state crosses over smoothly from the (3,3,1) state 
to the spin-polarized Pfaffian (overlap of < 0.92). How- 
ever, the increase in Ag^s, somewhat countcrintuitivcly, 
does not immediately destroy the Halperin state, but at 
first even increases the overlap. 

Finally, Fig. 5] shows the results of our Monte-Carlo 
study of the energies of the (3,3,1) and Pfaffian states. 
The correlation energies of both states were obtained 
from the finite size scaling of systems with N = 6 to 
N = 18 electrons as described above in Section IIV CI 
All data was obtained in Montc-Carlo simulations with 
10 7 samples. The uncertainty in the energy of the two- 
component states was obtained as the difference between 
linear and quadratic extrapolation of the background en- 
ergies (Fig. 2]), as this was larger than the bare numer- 
ical errors of the simulation. The energetic competi- 
tion of these two phases qualitatively recovers the pic- 
ture gained from studying the overlaps with the exact 
ground-state. Again, some finite amount of tunneling is 
required for the single component paired state to outcom- 
pete the Halperin state. As shown in Fig. [H the critical 
tunneling value A C SAS above which the Pfaffian state is 
energetically favored has a similar upturning shape as the 
boundary of large overlaps for the Pfaffian state in Fig. [3] 
However, there are some quantitative differences at small 
w, where the thermodynamic values indicate that polar- 
ization occurs at smaller values of the tunneling. 
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FIG. 5: (color online) Overlap between the exact Coulomb 
state of the quantum well for N = 6 particles at v = 1/4 
with the Halperin (5,5,3) state (a), the Pfaffian (b) and HR 
(c) states. The expectation value of the S z component of the 
pseudospin is plotted in (d). 



FIG. 6: (color online) Overlap between the exact Coulomb 
state of the quantum well for N = 8 particles at v = 1/4 
with the Halperin (5,5,3) state (a), the Pfaffian (b) and HR 
(c) states. The expectation value of the S z component of the 
pseudospin is plotted in (d). 



B. v = 1/4 in a quantum well 

We proceed with analyzing the quantum well at v = 
1/4 (Figs. [5]-[8|). Because of the rapid increase in size of 
the Hilbert space, there are only two system sizes acces- 
sible in ED at this filling factor: N = 6 and N = 8. The 
dimension of the L z = sector of the Hilbert space of the 
latter, taking into account discrete L z — > — L z symme- 
try, is on the order of 13 million, thus making N = 6 the 
only case amenable to study in great detail. However, for 
N = 6 we also must keep in mind the aliasing problem 
that occurs for (5,5,3) and the Pfaffian (there is no such 
problem for the HR state). We will present results for 
both particle numbers because of the important differ- 
ences between them. In view of the comments in section 
IIII1 wc note that the overlap with the (7,7,1) state is 
negligible in the range of widths w/Ib ^ 10.0, and there- 
fore we will exclude it from the present discussion of ED 
results. Note that, similarly to the (3,3,1) state in Sec. 
IV A| the (5,5,3) state hereinafter is defined as an eigen- 
state of the S x operator (if defined as an eigenstate of S z , 
the overlap with the exact ground state is negligible). 

In Fig. [S]we plot the overlap between the ground state 
of the quantum well for TV = 6 particles at v = 1/4 and 
the Halperin (5,5,3) state (a), the Pfaffian (b) and the 
HR state (c), accompanied by the expectation value of 
the S z component of the pseudospin (d). These results 



are reminiscent of v = 1/2 (Fig. [3]); however, due to 
the smaller energy scale and the gap, it is much easier 
to polarize the system at v = 1/4. For intermediate 
values of the width and small tunneling, the maximum 
overlap with the Halperin (5, 5, 3) state is high (0.96), but 
the region that would correspond to this phase is quite 
narrow in comparison to that of (3,3,1). On the other 
hand, the Pfaffian phase is much more extended. Given 
the intrinsic tunneling^ of the samples, which is of the 
order of Asas /(e 2 /ds) < 0.1, it seems more likely that 
the system will be found in this phase than the (5, 5, 3). 

The small island where the overlap abruptly goes to 
zero for large w/Ib is due to the ground state belonging 
to a sector with L > — this can be due to the admixture 
of compressible physics at large widths. The HR state ap- 
pears to be present in the transition region between one- 
and two-component phases, its overlap steadily increas- 
ing with w and peaking at 0.7 for w/Ib = 6.0. Because 
of the fact that the HR state occurs at a different shift 
on the sphere, we stress that the overlap presented here 
does not constitute a proof that it is an intermediary 
phase (moreover, the overlap drops rapidly when larger 
systems arc considered, sec Fig. [7J|. 

In Fig. [S] we plot the same quantities for the sys- 
tem of N = 8 particles which is expected to display 
weaker finite-size effects and docs not suffer from the 
aliasing problem. The (5,5,3) state is found in a siz- 
able parameter range, but the maximum overlap is mod- 
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FIG. 7: (color online) Overlap between the exact Coulomb 
state of the quantum well for N = 8 particles at v = 1/4 
and w/Ib = 10.5 with the Halperin (5,5,3) state, the Pfaffian 
and HR states (left axis). The expectation value of the S z 
component of the pseudospin is given on the right axis. The 
shaded region denotes where the ground state is no longer 
rotationally invariant (L > 0). 



erate compared to the case previously studied (0.74 for 
w/Ib — 4.5). While the HR state generally has a small 
overlap (not exceeding 0.2) and the evolution of (S z ) re- 
mains smooth, the striking difference in comparison with 
the N = 6 results (Fig. [5]) is the Pfaffian phase. Although 
it similarly develops with the increase in Ag^s, once the 
system reaches full polarization, the phase is destroyed. 

To shed more light on how this occurs, it is useful to 
look at the "cross section" of Fig.[B]for a fixed value of the 
width w/Ib = 10.5, chosen to represent the region where 
the Pfaffian phase is most clearly pronounced (Fig. [7]). 
Although the Pfaffian overlap peaks in the region where 
(5, 5, 3) starts to drop, very abruptly both overlaps fall to 
zero, and the ground state is no longer rotationally invari- 
ant. The fact that L > is a hallmark of compressibility. 
Precisely at the transition point, a small kink is now vis- 
ible in (S z ). The origin of this kink or the reason why 
the ground state is obtained in L > sector are not en- 
tirely clear at present. However, the zero overlap with the 
Pfaffian beyond As as / '(e 2 / 'tis) = 0.1 (where the ground 
state reduces to a spinless case) agrees with our results of 
the section HH Notice that a compressible ground state 
with L > may also indicate a phase with modulated 
charge density, such as the Wigner crystal. Indeed, an 
insulating behavior, as one would expect for an electron 
crystal, has been found at filling factors slightly above 
v = 1/541 Such a state is not captured in the present 
ED calculations on the sphere, and the question whether 
a Wigner crystal is the true ground state at large values 
of A.sas in a wide quantum well is beyond the scope of 
the present paper. 

We refer to Monte-Carlo simulations (Fig. [5]) to ob- 
tain additional information about the candidate incom- 
pressible states from larger model systems. We include 
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FIG. 8: (color online) Results from our Monte-Carlo study 
of states at v = 1/4: (a) correlation energies of the 
Pfaffian-, (5,5,3)- and (7,7,l)-states with respect to the single- 
component background, as a function of the well-width w in 
the thermodynamic limit; (b) difference in energy between 
the different Halperin states and, in particular, the tunneling 
strength Ag AiS for the Pfaffian state to be favoured over the 
(5,5,3)-state, and (c) correlation energies for the Pfaffian state 
in units of rescaled magnetic length^ l' B = \N / {vN^Y^Ib 
for some values of w. note the particularly high value at 
TV = 8. 



systems with N = 6 to N = 16 electrons in the finite 
size scaling for the groundstate energies, again using 10 7 
Monte-Carlo samples, and taking errors as the difference 
between linear and quadratic extrapolation of the back- 
ground energies. The results of this study are summa- 
rized in Fig. where we compare the Pfaffian to the 
(5,5,3) and (7,7,1) Halperin wave functions. Again, a 
two-component state is always preferred in the absence of 
tunneling. At the layer separations shown, this is (5,5,3) 
as shown in Fig. [8ja). These data also confirm that the 
(7,7,1) state becomes relevant only at large well width 
W > 101b- In Fig. [Hb), we display the value of tunnel- 
ing A C SAS required to polarize the system into the paired 
Pfaffian phase. This feature of the energetic competi- 
tion of (5,5,3) and the Pfaffian is very close to the results 
obtained in ED for N = 6 both qualitatively and quan- 
titatively: the shape of A c SAS (w) is nearly linear and 
reproduces the location where the overlaps with the ex- 
act ground-state cross over between the two trial states, 
as was shown in Fig. [5l The splitting A$as required for 
the Pfaffian to be the ground state is significant, and 
probably larger than the splitting in the experiments 
of Luhman et alr^, which can be estimated to about 
Asas ~ 0.069e 2 /ds at the sample width w as 101b and 
their baseline electron density. 

This similarity between the energetics in the thermo- 
dynamic limit and the exact spectrum for N = 6 particles 
may be circumstantial. However, there is another indi- 
cation that the very different behavior at N = 8 might 
be exceptional. In Fig. EJc), we show the correlation en- 
ergies of the Pfaffian state for different system sizes N 
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and well widths w. This representation reveals the case 
of TV = 8 as having particularly high energy. This may 
be a finite-size effect that can be explained in the com- 
posite fermion picture. The Pfaffian wavefunction can 
be expressed as a paired state of 4 CF feeling one quan- 
tum of negative effective flux . 22 ' 42 The shell structure of 
these composite fcrmions on the sphere yields filled shell 
states for TV = 6 and TV = 12, whereas for TV = 8 two 
CFs remain in the highest, partially filled shell. In this 
configuration, CFs are susceptible to follow Hund's rule 
by maximising their angular momentum, breaking rota- 
tional invariancc. 

For TV = 8 and TV = 10, Hund's rule predicts an an- 
gular momentum of L = 4, which is indeed found in ED. 
This gives us confidence that the system is still described 
by liquid-like composite fermion physics at large As^s- 
We therefore consider the competition between a Hund's 
rule state and the paired Pfaffian state. For a similar 
situation with weak pairing in a v = 1/2 + 1/2 bilayer 
system at large layer separation, it was argued^ that for 
larger systems, the shell-filling effects, and Hund's rule, 
should become less important whereas the pairing effects 
will remain the same strength, as only ~ VTV benefit from 
Hund's rule, whereas all (~ TV) particles within some gap 
energy of the Fermi surface contribute to pairing. 

Although the above argument speaks in favor of the 
possibility for a paired Pfaffian state to be realized at 
v = 1/4 for large tunneling gap Ag^s, we insist on 
the variational character of the Monte-Carlo calculations. 
In these calculations, we have indeed considered several 
competing candidate wave functions for a liquid ground 
state at this filling factor. However, this analysis may not 
eliminate the possibility that a compressible state, such 
as that seen in ED, or even other incompressible phases 
may indeed be singled out as a true ground state of the 
system. 

Finally, we would like to point out that in ED it is 
possible to calculate the quantity that we refer to as the 
"charge gap" , 

AE = Eif^+i + EV.a^-i — 2£V jA r #) (18) 

where En,n^ is the ground state energy for a given num- 
ber of particles TV and number of flux quanta TV^. This 
quantity probes the response of the system to the intro- 
duction of quasiparticles / quasiholes on top of the ground 
state, and its dependence on Ag^s has been used to de- 
lineate between the one-component and two-component 
phases^ With the appropriate finite-size corrections, 
Eq. (|T8j) should correspond to the experimentally mea- 
surable "activation" gap 1 ^ that governs the temperature 
scaling of longitudinal resistance R xx ~ exp (— AE/2T). 
For states that undergo a typical one- to two-component 
transition, such as the one at v = 2/3 (for small tun- 
neling, it is the state of two decoupled Laughlin liquids, 
f = 1/3 + 1/3, which develops into a single-component 
2/3-state for large tunneling amplitudes 4 ^), the charge 
gap (|18|) displays a minimum as a function of As as in 
the center of the transition regionJ^ On the other hand, 



for v = 1/2 where the tunneling-driven transition con- 
nects the (3,3, 1) state and the Pfaffian, the charge gap 
(fT5|) shows an upward cusp. Our calculations of the 
charge gap Q18p in the case of v = 1/4 indicate that this 
quantity is a less robust way to characterize the nature 
of the ground state than the calculation of the overlaps 
with trial wave functions. While for TV = 6 particles at 
v = 1/4 the charge gap displays a minimum as a func- 
tion of Asas, there is a very weak dependence of AE on 
As as when a larger system of TV = 8 particles is consid- 
ered. Thus finite-size effects are too strong in order to 
extract useful information from Eq. (|18[) in small systems 
that can be treated by ED. 



VI. CONCLUSION 

In this paper we have presented a systematic study of 
several candidates for the ground-state wave function at 
the recently observed 1 ^ fraction v = 1/4. Assuming that 
the (pseudo)spin plays no role, i.e., in a one-component 
picture, the generalized Moore-Read Pfaffian state ([1]) 
shows high overlap for the values of the sample width 
which are on the order of those in the experiment of 
Ref. [H, but only if the confinement in the perpendic- 
ular direction is modeled by the ZDS model (J3J). For 
other confinement models (J2HU) it was not possible to re- 
produce such high values of the overlap. We believe that 
this inconsistency means that the high overlap must be 
due to a special softening of the pseudopotentials that 
occurs as a pathology of ZDS model but does not appear 
in other (more realistic) confinement models. 

Therefore, the existence of a fractional quantum Hall 
state at v = 1/4 is necessarily linked to the specific fea- 
tures of the quantum well used in Ref. [TH that enable 
the multicomponent physics to manifest itself. Addi- 
tional degrees of freedom in our theoretical study are 
conveniently taken into account within the quantum-well 
model, which is the simplest model that can naturally 
interpolate between a single layer and bilayer charge dis- 
tribution as the parameters w and As^s are varied. This 
two-parameter model is related to the previous studies 1 ^ 
of the true bilayer with tunneling at v = 1/2 (which 
had to assume at least three independent parameters) by 
reproducing the same physical picture of the crossover 
between the (3,3,1) state and the Pfaffian. 

At the filling factor v = 1/4, we have not been able 
to produce clear cut evidence for the expected crossover 
between the (5,5,3) state and the Pfaffian in ED, due to 
the strong finite-size effects in case of the latter. We have 
shown that the (5,5,3) state is indeed present for a range 
of widths and small tunneling gaps As^s, but its max- 
imum overlap is not as high as that of the (3,3,1) state. 
ED cannot delimit the range of parameters for the Pfaf- 
fian phase due to the difference in the results for the two 
available system sizes, TV = 6 and TV = 8, and the effect 
of compressible physics which is difficult to treat within 
the spherical geometry. However, our Monte-Carlo simu- 



12 



lations go some way towards clarifying the situation. The 
correlation energies of the Pfaffian state reveal N = 8 as 
a particularly unfavorable system size. We can explain 
this from the finite-size effect in terms of filling shells of 
CF orbitals on the sphere. The competing L ^ states 
at N — 8, as well as N = 10, seem to be related to Hund's 
rule for CFs. However, the competition between Hund's 
rule and pairing is likely favorable for the paired state in 
the thermodynamic limit. In addition, projecting from 
the two-component model onto the fully polarized (spin- 
less) case, on the other hand, can be seen as analogous to 
the scenario of LL mixing^ which may provide another 
mechanism to stabilize the Pfaffian state via generating 
three-body terms in the effective interaction. Such effects 
are beyond the scope of present paper. By analyzing the 
competition between the paired single component and 
the Halpcrin states from their variational wave functions, 
we find, in the Monte-Carlo simulations, that the tunnel- 
ing gap A C SAS required to form a single-component state 
roughly behaves linearly as 1.0 x (w/Ib) x 10 -2 e 2 /e?B- 
Although the tunneling splitting indicated for the exper- 
iment described in Ref. Il6| is not far below the transi- 
tion between the Pfaffian and (5,5,3), our numerics still 
show it safely in the two-component regime of the (5,5,3)- 
wavefunction. 



Although we believe that our quantum well model 
takes properly into account the effects of finite thick- 
ness, we have entirely neglected the effect of the in- 
plane magnetic field which may nevertheless prove es- 
sential in order to stabilize the incompressible state at 
v = 1/4. The existing experimental work^S on the 
v — 2/3 state witnessed that the introduction of an in- 
plane magnetic field may lead to a strengthening of the 
minimum in R xx , thus inducing the same one-component 
to two-component transition as by varying Ag^g. Simi- 
lar strengthening occurs for v = 1/2 if the tilt is not too 
large.— Therefore, the application of the in-plane field 
may be a likely reason to further stabilize the (5,5,3) state 
at v = 1/4 if the symmetric-antisymmetric gap As^g is 
sufficiently small. However, Ref. [ly also pointed out the 
difference between v = 1/2 state and v = 1/4 state: 
when the electron density is increased, the former dis- 
plays a deeper minimum in R xx while the latter remains 
largely unaffected. This difference suggests that in the 
case of v = 1/4 the quantum- well ground state may be 
effectively fully polarized and in the class of the Pfaffian 
rather than the two-component, (5,5,3) state. 

In order to answer without ambiguity which of the two 
possibilities is actually realized in the quantum well un- 
der the experimental conditions of Ref. [l6l it would be 
useful to know the dependence of the activation gap as 
a function of Ag^s and also as a function of transferred 
charge from the front to the back of the quantum well us- 
ing a gate biasing. These results would help to discrimi- 
nate between the one-component and two-component na- 
ture of the ground state. 
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APPENDIX A: EFFECTIVE BILAYER 
DESCRIPTION OF THE WIDE QUANTUM 
WELL 



As in Sec. IIV1 wc consider the quantum well to be 
symmetric around w/2, i.e. the lowest subband (f) state 
is symmetric, and the first excited one (J.) is antisym- 
metric. Furthermore, we consider, in this section, the 
electrons to be in the 2D plane, for illustration reasons, 
although the conclusions remain valid also in the spher- 
ical geometry. In this appendix, wc yield the derivation 
of the effective bilaycr description of the wide quantum 
well. 

The interaction part of the Hamiltonian §§§ consists 
of a density-density interaction and terms beyond, which 
may be described as a spin-spin interaction. Indeed, the 
density-density part consists of the effective interactions 

Vi m . Notice that the 



I V^ n , V^ Ll , and V^ U 



interactions in the first excited subband (J,) are generally 
weaker than in the lowest one (J.) because the wave func- 
tion §5§ <pi(z) possesses a node at w/2, in the center of 
the well, i.e. > V^ U . With the help of the (spin) 

density operators (|11|) and (TT2")) . the density-density part 
of the interaction Hamiltonian reads 



H 



2 v su(2) (q)p(-qMq) 

q 

+2j2V s U^S z (~q)S z (q) 
q 

+ 5>i(q)p(-q)S*(q) , 
q 

in terms of the SU(2)-symmetric interaction 



(Al) 



l r 

4 



Vsu&M = 7 ^D TTT (q) + ^D u (q) + 2^ 2 T D m (q) 



(A2) 



and the SU(2)-symmetry breaking interaction terms 



VM = 7 ^2D TT (q) + V^ U (q) - 2V^ U (q) (A3) 



and 



VS (q) 



^D m (q)-^D (q) 



(A4) 
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The remaining 12 interaction terms, which may not 
be treated as density-density interactions, fall into two 
different classes; the 8 terms with three equal spin ori- 
entations a and one opposite —a are zero due to the 
antisymmetry of the integrand in Eq. (|10[) . The remain- 
ing 4 interaction terms with two |-spins and two |-spins 
are all equal due to the symmetry of the quantum well 
around w/2, 

\r x — T/TTil _ T/iiTT _ T/T11T _ T/lTTi / a c\ 

v sb = V 2D ~ V 2D — V 2D ~ V 2D • \ A0 ) 

They yield the term 

H' th = 2 £ ^(q)^(-q)^(q) , (A6) 
q 

I 



which needs to be added to the interaction Hamiltonian 
(|Alj) . as well as the term 



H SAS = -A S AsS z (q = 0), (A7) 



which accounts for the electronic subband gap between 
the t and the J, levels. 

Collecting all terms, the Hamiltonian ([9]) thus becomes 



H = 2E V)(#H^ + 2^](q)?(-q)S l (q) + 2^^(q)^(-q)S'(q) 
q q q 

+ ^VJ(q)p(-q)5 2 (q) - A SAS 5 z (q - 0) . (A8) 



Several comments are to be made with respect to this re- 
sult. First, we have checked that for the infinite-square- 
well model as well as for a model with a parabolic confine- 
ment potential there is a natural hierarchy of the energy 
scales in the Hamiltonian (|A8|) . 

V SU (2) > V% > VJ > VA . (A9) 

This hierarchy is valid both for the interaction potentials 
in Fourier space as for the pseudopotentials. 

Whereas the first term of the Hamiltonian describes 
the SU(2)-symmetric interaction, the second and the 
third one break this SU(2) symmetry. Because V^(q) > 
V* b (q) > for all values of q, states with no polarization 
in the x- and the ^-direction are favored, with (S x ) = 
and (S z ) = 0, respectively. Due to the hierarchy (|A9|) of 
energy scales, a depolarization in the x-direction is more 
relevant than that in the z-dircction. These terms are 
similar to those one encounters in the case of a bilayer 
quantum Hall system, where due to the finite layer sepa- 
ration a polarization of the layer isospin in the z-direction 
costs capacitive energy^ 



The fourth term of the Hamiltonian (|A8[) is due to the 
stronger electron-electron repulsion in the lowest elec- 
tronic subband as compared to the first excited one, 
where the wave function possesses a node at z = w/2. 
In order to visualize its effect, one may treat the den- 
sity, which we consider to be homogeneous in an incom- 
pressible state, on the mean-field level, (p(q)) = ^<5q,o, 
in which case the fourth term of Eq. (|A8|) becomes 
^^(q = 0)<^ z (q = 0) and, thus, has the same form 
as the subband-gap term (|A7j) . It therefore rcnormalizcs 
the energy gap between the lowest and the first excited 
electronic subbands, and it is natural to define the effec- 
tive subband gap as 

A s A s -> A s as = A«jas - ^Vg(q = 0) 

9 

6 W 

= Asas-Iv-tt , (A10) 
els IB 

where 7 is a numerical prcfactor that depends on the 
precise nature of the well model. 



1 W. Pan, J. S. Xia, H. L. Stormer, D. C. Tsui, C. Vicente, 
E. D. Adams, N. S. Sullivan, L. N. Pfeiffer, K. W. Baldwin, 
and K. W. West, Phys. Rev. B 77, 075307 (2008). 

2 R. B. Laughlin, Phys. Rev. Lett. 50, 1395 (1983). 

3 J. Jain, Composite Fermions (Cambridge University Press, 
Cambridge, England, 2007). 

4 F. D. M. Haldane, Phys. Rev. Lett. 51, 605 (1983). 

5 R. Willet, Phys. Rev. Lett. 59, 1776 (1987). 



6 F. D. M. Haldane and E. Rezayi, Phys. Rev. Lett. 60, 956 
(1988). 

7 G. Moore and N. Read, Nucl.Phys. B 360, 362 (1991). 

8 N. Read and E. Rezayi, Phys. Rev. B 54, 16864 (1996). 

9 E. H. Rezayi and F. D. M. Haldane, Bull. Am. Phys. Soc. 
32, 892 (1987). 

J. Eisenstein in Perspectives in Quantum Hall Effects, 
edited by S. Das Sarma and A. Pinczuk (John Wiley & 



14 



Sons, New York, 1997). 

11 Y. W. Suen, L. W. Engel, M. B. Santos, M. Shayegan, and 
D. C. Tsui, Phys. Rev. Lett. 68, 1379 (1992). 

12 Y. W. Suen, H. C. Manoharan, X. Ying, M. B. Santos, and 
M. Shayegan, Phys. Rev. Lett. 72, 3405 (1994). 

13 S. He, S. Das Sarma, and X. C. Xie, Phys. Rev. B 47, 4394 
(1993). 

14 K. Nomura and D.Yoshioka, J. Phys. Soc. Jpn. 73, 2612 
(2004). 

15 B. I. Halperin, Helv. Phys. Acta 56, 75 (1983). 

16 D. R. Luhman, W. Pan, D. C. Tsui, L. N. Pfeiffer, K. W. 
Baldwin, and K. W. West, Phys. Rev. Lett. 101, 266804 
(2008). 

17 F. D. M. Haldane and E. H. Rezayi, Phys. Rev. Lett. 54, 
237 (1985). 

18 B. Andrei Bernevig and F. D. M. Haldane, Phys. Rev. Lett. 

100, 246802 (2008). 

19 F. C. Zhang and S. Das Sarma, Phys. Rev. B 33, 2903 
(1986). 

20 M. Peterson, Th. Jolicoeur and S. Das Sarma, Phys. Rev. 
B 78, 155308 (2008). 

21 K. Park, V. Melik-Alaverdian, N. E. Bonesteel, and J. K. 
Jain, Phys. Rev. B 58, R10167 (1998). 

22 G. Moller and S. H. Simon, Phys. Rev. B 77, 075319 
(2008). 

23 E. Rezayi and F. D. M. Haldane, Phys. Rev. Lett. 84, 4685 
(2000). 

24 Z. Papic and M. V. Milovanovic, Phys. Rev. B 75, 195304 
(2007). 

25 G. Moller, S. H. Simon and E. H. Rezayi, Phys. Rev. Lett. 

101, 176803 (2008). 

26 G. Moller, S. H. Simon and E. H. Rezayi, Phys. Rev. B 
79, 125106 (2009). 

27 M. O. Goerbig and N. Regnault, Phys. Rev. B 75, 241405 
(2007). 

28 R. de Gail, N. Regnault, and M. O. Goerbig, Phys. Rev. 
B 77, 165310 (2008). 

29 N. Regnault, M. O. Goerbig and Th. Jolicoeur, Phys. Rev. 
Lett. 101, 066803 (2008). 

30 D. Yoshioka, A. H. MacDonald, and S. M. Girvin, Phys. 
Rev. B 39, 1932 (1989). 



31 H. C. Manoharan, Y. W. Suen, T. S. Lay, M. B. Santos, 
and M. Shayegan, Phys. Rev. Lett. 79, 2722 (1997). 

32 E. H. Rezayi and N. Read, Phys. Rev. Lett. 72, 900 (1994). 

33 M. V. Milovanovic and N. Read, Phys. Rev. B 53, 13559 
(1996). 

34 M. V. Milovanovic, Th. Jolicoeur and I. Vidanovic, 
larXiv:0902.1719l 

35 M. Abolfath, L. Belkhir, and N. Nafari, Phys. Rev. B 55, 
10643 (1997). 

36 For a review, see K. Moon, H. Mori, K. Yang, S. M. Girvin, 
A. H. MacDonald, I. Zheng, D. Yoshioka and S.-C. Zhang, 
Phys. Rev. B 51, 5138 (1995); S. M. Girvin and A. H. Mac- 
Donald, in Perspectives in Quantum Hall Effects, edited by 
S. Das Sarma and A. Pinczuk (John Wiley & Sons, New 
York, 1997). 

37 C. Toke and J. K. Jain, Phys. Rev. Lett. 96, 246805 (2006). 

38 G. Fano, F. Ortolani and E. Colombo, Phys. Rev. B 34, 
2670 (1986). 

39 R. Morf, N. d'Ambrumenil, and B. I. Halperin, Phys. Rev. 
B 34, 3037 (1986). 

40 T. S. Lay, T. Jungwirth, L. Smrcka, and M. Shayegan, 
Phys. Rev. B 56, R7092 (1997). 

41 H. Jiang, R. L. Willet, H. L. Stormer, D. C. Tsui, L. N. 
Pfeiffer et K. W. West, Phys. Rev. Lett. 65, 633 (1990); 
H. Jiang, H. L. Stormer, D. C. Tsui, L. N. Pfeiffer et K. 
W. West, Phys. Rev. B 44, 8107 (1991). 

42 G. Moller and S. H. Simon, Phys. Rev. B 72, 045344 
(2005). 

43 C. T6ke, N. Regnault and J. K. Jain, Sol. Stat. Com. 144 
(12) 504 (2007). 

44 In practice, the use the effective potential (|17[1 is limited 
to moderately large systems with N$ < 60. For larger sys- 
tems, it is more suitable to use effective potentials based 
explicitly on the asymptotic behaviour of the in-plane in- 
teraction i|iop. 

45 A single component v = 2/3 state can be regarded either as 
the particle- hole conjugate of the v = 1/3 Laughlin state, 
or as composite fermions at negative effective flux filling 
p = -2 CF-LLs.i 2 - 



